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The paramagnetic properties in non-centrosymmetric superconductors with and without an- 
tiferromagnetic (AFM) order are investigated with focus on the heavy Fermion superconduc- 
tors, CePtsSi, CeRhSia and CelrSis. First, we investigate the spin susceptibility in the linear 
response regime and elucidate the role of AFM order. The spin susceptibility at T = is in- 
dependent of the pairing symmetry and increases in the AFM state. Second, the non-linear 
response to the magnetic field are investigated on the basis of an effective model for CePtsSi 
which may be also applicable to CeRhSia and CelrSia. The role of antisymmetric spin-orbit 
coupling (ASOC), helical superconductivity, anisotropic Fermi surfaces and AFM order are 
examined in the dominantly s-, p- and d-wave states. We emphasize the qualitatively impor- 
tant role of the mixing of superconducting (SC) order parameters in the p-wave state which 
enhances the spin susceptibility and suppresses paramagnetic depairing effect in a significant 
way. Therefore, the dominantly p-wave superconductivity admixed with the s-wave order pa- 
rameter is consistent with the paramagnetic properties of CePtaSi at ambient pressure. We 
propose some experiments which can elucidate the novel pairing states in CePtaSi as well as 
CeRhSia and CelrSia. 

KEYWORDS: Superconductivity without inversion center; antiferromagnetic superconductor; Pauli para- 
magnetic effect; spin susceptibility 



1. Introduction 

Since the discovery of superconductivity in the non- 
centrosymmetric heavy Fermion compound CePtsSi,^'^ 
superconductivity in materials without inversion cen- 
ter is attracting growing interest. Many new non- 
centrosymmetric superconductors (NCSC) with un- 
usual properties have been identified among heavy 
fermion systems such as Ulr,"^ CeRhSis,^'^ CelrSia,^''' 
CeCoGes^ and others like Li2Pd:,Pt3_^B,9 YaC 



2^3, 



RhzGag, IrzGag,"'!^ MgioIrigBig," ReaW^'' and some 
organic materials. The aspects of missing inversion 
symmetry are also of great interest for other materials. 
For example, the spin Hall effect in the semiconductor^^ 
and the helical magnetism in MnSi^^ are very active re- 
search fields. 

NCSC adds several unusual aspects to the properties of 
superconductivity. One immediate consequence of non- 
centrosymmetricity is the necessity for an extended clas- 
sification scheme of Cooper pairing states, as parity is not 
available as a distinguishing symmetry. Using the tradi- 
tional scheme the SC states here may be represented as 
a mixture of pairing states of even and odd parity, or, 
equivalently, their spin configuration is a superposition 
of a singlet and a triplet component. This is a conse- 
quence of the presence of antisymmetric spin-orbit cou- 
pling (ASOC) in non-centrosymmetric materials. Re- 
cent theoretical studies led to the discussion of various in- 
triguing properties which could appear in NCSC, such as 
the magneto-electric effect, ""^^"^^ the unusual anisotropic 
spin susceptibility, 22-28 ^j-^g occurrence of an anoma- 
lous coherence effect in NMR l/TiT,2i'29 ^jj^g unusual 



origin of nodes in the SC gap, 2'^' 29-32 -j-j^g realization of 
the helical SC phase, the possible appearance of 
Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state at zero 
magnetic field, '^^ de Haas-van Alphen effect, *° various 
novel impurity effects, ^^"''^ and vortex core states^^ and 
unconventional features in quasiparticle tunneling and 
Josephson effect. ^^"^^ 

The non-centrosymmetric heavy fermion superconduc- 
tors, e.g. CePtaSi, Ulr, CeRhSis, CelrSis and CeCoGes 
are of particular interests because the Cooper pairing 
is most likely unconventional (non-s-wave) due to the 
strong electron correlation. Although many studies have 
been devoted to this topics, there is no consensus on 
the symmetry of pairing in these compounds so far. The 
symmetry of Cooper pairs may be determined by the 
paramagnetic properties such as the spin susceptibility 
below Tc- 

In centrosymmetric superconductor, the spin suscep- 
tibility is a distinguishing feature for the spin config- 
uration of the pairing state, as it decreases below Tc 
for the spin singlet superconductor and remains con- 
stant in the case of spin triplet pairing, if the mag- 
netic field is perpendicular to the d-vector (parallel to 
the equal-spin direction). The measurements of the 
Knight shift which is proportional to the spin suscep- 
tibility have played an important role for the identifi- 
cation of SC state in various compounds. For super- 
conductors with very high Hc2 probing effects of para- 
magnetic limiting can give also insight into the pair- 
ing symmetry and has been applied in connection with 
NCSC. However, the response to the magnetic field is 
not so straightforward in non-centrosymmetric systems. 
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As mentioned above, spin singlet and triplet compo- 
nents are mixed in the pairing state. Furthermore the 
band splitting induced by the ASOC affects the mag- 
netic properties. Therefore, it is necessary to clarify 
the magnetic properties very carefully before drawing 
strong conclusions. In this context also the influence 
of AFM order on the (magnetic) properties of the SC 
phase is an important point to investigate, since all 
presently known non-centrosymmetric heavy Fermion su- 
perconductors, i.e. CcPtaSi, Ulr, CeRhSia, CelrSia and 
CeCoGes, coexist with the magnetism. In CePtsSi at 
ambient pressure, superconductivity (Tc = 0.75K) co- 
exists with AFM order (Tn = 2.2K).i'5^'" The AFM 
order can be suppressed by pressure and vanishing at 
the critical value of P ~ 0.6GPa. The SC phase is more 
robust and a purely SC phase exists beyond the criti- 
cal pressure, (P > 0.6GPa).58-60 CeRhSis,"'^ CelrSig'''^ 
and CeCoGes^ are AFM at ambient pressure and super- 
conductivity appears only under substantial pressure. Al- 
though most of the theoretical studies except for Refs. 27, 
30 and 61 neglected the AFM order so far, it turns 
out that the magnetism affects the electronic state pro- 
foundly. It has been shown that a gap line-node behavior 
could be induced by the AFM order for the pairing state 
with (lominantly p-wavc componcnt,^^'"'"'^^ which may 
explain the experimental results in CePtsSi at ambient 
pressure.^^"^^ 

In this paper we investigate the linear as well as the 
non-linear response regime of the NCSC in a magnetic 
field with the aim to provide guidelines to identify the 
pairing symmetry based on magnetic properties. Before 
going into details we briefly summarize the main conclu- 
sions of our study. It is known that in the linear response 
regime the paramagnetic properties are universal i.e. the 
spin susceptibility is independent of the pairing symme- 
try. In the presence of Rashba-type ASOC, the spin sus- 
ceptibility along the c-axis is constant through Tc while it 
decreases along the a6-plane to half of the normal state 
value at T = 0, in absence of AFM order. 22-28 The 
influence of lielicdty (Cooper pairs possess a finite mo- 
mentum) in NCSC on the behavior of the susceptibility 
turns out to be negligible. On the other hand, the fold- 
ing of Brillouin zone due to the AFM order significantly 
affects the spin susceptibility in the SC phase. The spin 
susceptibility for the magnetic field perpendicular (par- 
allel) to the staggered moment is increased (decreased) 
by the AFM order.^^ 

The non-linear response to the magnetic field is im- 
portant when the magnetic field is comparable to or 
higher than the standard paramagnetic limiting field 
Hp ~ 1.2kBTc/nB- It should be noted that most of the 
experimental studies, such as the Knight shift and criti- 
cal magnetic field i?c2, have been carried out in the non- 
linear response region. ^'2'^-'''^8, 66-70 rpj^g pairing state in 

NCSC can be identified by the measurements in the non- 
linear response regime because the paramagnetic prop- 
erties depend on the pairing symmetry in contrast to the 
situation in the linear response regime. 

We show that the critical magnetic field ifc2 along the 
a6-plane is significantly enhanced in the non-linear re- 
sponse regime by the formation of helical SC state. This 



enhancement coincides with the non-linear increase of 
the helicity of the SC order parameter. Hc2 furthermore 
rises for the dominantly p-wave state owing to the mix- 
ing of SC order parameters. These effects, namely (i) the 
formation of the helical SC state and (ii) the mixing of 
SC order parameters, are quantitatively important for 
anisotropic Fermi surfaces. AFM order significantly en- 
hances the effect (ii) and also boosts Hc2- In this case, 
the spin susceptibility remains nearly constant through 
Tc. On the other hand, these effects are negligible in the 
dominantly spin singlet pairing state. Since the infiuence 
of AFM order is quantitatively important, the paramag- 
netic properties of the SC phase in the AFM state pro- 
vide a means to distinguish between pairing states with 
dominant spin triplet and singlet component. 

Among the non-centrosymmetric heavy fermion super- 
conductors, CePtsSi has been investigated in most detail 
because the superconductivity exists at ambient pressure 
while others require substantial pressure to become su- 
perconducting. Therefore, we pay particular attention to 
the situation in CePtaSi, and discuss the pairing sym- 
metry by comparing the experiments'"^' 2' '^■^^'^^'^^ with 
our theoretical results. The paramagnetic properties of 
CePtaSi look puzzling at first sight because the experi- 
mental results are incompatible with the theoretical re- 
sults within the linear response theory and without tak- 
ing into account the AFM order.^^' 22-26, 28 jj^ q^j. present 

study we show that the experimental results are consis- 
tent with the theoretical results for the dominantly p- 
wave state by taking into account the AFM order as well 
as the non-linear response to the magnetic field. 

Moreover we propose further test experiments which 
could strengthen our conclusions. First, the influence of 
AFM order can be examined by the pressure which sup- 
presses the AFM order. Second, the 2-fold anisotropy 
in the a6-plane arises from the AFM order and the 
anisotropy is qualitatively different between the (lomi- 
nantly p-wave, inter-plane d-wave and intra-plane s- or 
d-wave states. Future experimental studies of these kind 
could help identify the pairing symmetry in CePtaSi, 
CeRhSig, CelrSia and CeCoGes. 

The paper is organized as follows. In §2 we summarize 
the linear response theory for the paramagnetic proper- 
ties in NCSC. In §3 we introduce the effective model 
for CePtsSi which could be also applied to CeRhSia, 
CelrSia and CeCoGea. The paramagnetic properties in 
the magnetic field along the a6-plane are investigated in 
§4 and §5. The non-linear response to the magnetic field 
for the dominantly s-wave state is investigated in §4. In 
§5, which is the main part of this paper, we show the 
magnetic properties in the dominantly p-wave state. The 
influences of the helical superconductivity, anisotropic 
Fermi surface and AFM order are elucidated. The pairing 
symmetry of CePtaSi is discussed in §6 by comparing the 
experimental results with our theoretical results. Some 
test experiments are proposed for CePtaSi, CeRhSia and 
CelrSia in §7. In §8, nature of the helical SC state is 
investigated in details. We show the crossover from the 
helical SC state with long wave length to that with short 
wave length. These results are summarized and some dis- 
cussions are given in §9. 
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2. Linear Response Theory 

In this section wc investigate the linear response 
regime of the NCSC in a magnetic field, and study the 
magnetic properties in the paramagnetic (PM) and in the 
AFM phase. The latter wc consider both for the case of 
a centrosymmetric and a non-centrosymmetric system. 

2.1 General spin susceptibility 

In a first step we derive a general expression for the 
spin susceptibility in the SC state on the basis of the 
extended BCS Hamiltonian, given by 

H = H^ + Hso + Haf+Ha, (1) 

fe,5 

Hso = aJ29{k)-S{k), (3) 

k 

HAF = -Y,hQ-SQ{k), (4) 

k 

HA = -Y.}^^,s,Ak)cl,_,A,,s 
s,s' yk 

where k± = k ± qu/2, S{k) = Z^gg' ct^^ ct ^cj^, and 

Sqik) = J2ss' ^ss'ci_^^ ^cj:^^,. Here is the total mo- 
mentum of Cooper pairs. Note that gk is zero in the 
usual BCS state while that is finite in the helical SC 
state. ^^"'^^ In NCSC the helical SC state can be realized 
under magnetic field above Hd. We consider a tetrago- 
nal crystal lattice and assign the x-, y- and a^axis to Or, 
b- and c-axis, respectively. 

The first term in cq. (1) describes the dispersion rela- 
tion without ASOC and AFM order. In this subsection 
we do not identify the specific dispersion of the electrons 
and assume e(fc) as general. 

The second term Hgo describes the ASOC due to 
the lack of inversion symmetry. This term preserves 
time reversal symmetry, if the (^-vector is odd in fc, i.e. 
g{—k) — —g{k). We consider a Rashba-type spin-orbit 
coupling^^ as is realized in CePtsSi, CcRhSia, CcIrSis 
and CeCoGea.^'* Because the detailed momentum de- 
pendence of g{k) is unknown, we express it in terms of 
velocities v{k) = de{k)/dk : g{k) = (— Uy (fc), Wx(fc), 0)/w. 
This choice at least preserves the correct periodicity in k- 
space. The detailed form of the 5- vector is anyway unim- 
portant in the following. We normalize the (j-vector g{k) 
by the average velocity v [v"^ = ;^Z^fc^^x(^)^ + Vy(ky] 
so that the coupling constant a has the dimension of 
energy. Wc assume the relation |A,;.s_,,/(fc)| <§; |a| ^ Sp 
throughout this paper {sp is the Fermi energy). This re- 
lation is valid for the most of NCSC such as CePtsSi, 
Ulr, CeRhSig, CcIrSig and CcCoGcg. 

The third term Haf is taken into account to inves- 
tigate the role of AFM order which enters through the 
staggered field hq. We focus on A- type AFM order, i.e. 
ferromagnetic sheets in the a6-plane are staggered along 



the c-axis, giving rise to Q = (0, 0, tt). This spin structure 
is realized in GePtsSi^'' as well as the centrosymmetric 
superconductor UPd2Al3^^'''^ where the magnetic mo- 
ments are aligned in the a&-planc. A different AFM state 
has been reported for CeRhSis'^'' and the magnetic struc- 
ture is not clearly identified for CelrSis so far. However, 
the qualitative role of AFM order can be captured by 
Q = (0, 0, tt) in the simple cases. 

The last term describes the mean field term of the 
SC order. The order parameter is given by Ai^s,s'{k) and 
^2,3.8' (k). The second component A2.s,s' (k) only appears 
in the case of superconductivity coexisting with AFM 
order (A2,s,s'(^) = for hq = 0). The order parameter 
has both the spin singlet and triplet components owing 
to the ASOC. 

It is more transparent for the following discussion to 
consider the order parameter in the band basis because 
the superconductivity is mainly induced by the intra- 
band Cooper pairing when |A| ^ \a\. Ignoring the order 
parameters describing the inter-band pairing, we obtain 
the simplified Hamiltonian as, 
4 

7=1 k 

-[^7(^)<,_g_<g^+/i-c-], (6) 

where Xlfe' is restricted to the summation within j/czl < 
tt/2. The dispersion relation ej{k) takes into account the 
ASOC and AFM order and is obtained by the unitary 
transformation as, 

U\k)H{k)U{k) = {ei{k)5ij), (7) 
where the 4x4 matrix H{k) is expressed as, 

( .(^ S) ) ■ 

We define e{k) = e{k)a^^^ + ag{k)a and a represent 
the three Pauli matrices and (t^^^ is the 2 x 2-unit ma- 
trix. The four bands e-y{k) are non-degenerate except for 
the special momentum, if a 7^ 0. Moreover, the relation 
e^(— fc) = e-y{k) is hold owing to the time-reversal sym- 
metry. 

The order parameter is expressed in the band basis as, 
Aband(fc) = U\k+)^.Mk)U*{-k-), (9) 

where 

A ■ ik\-( ^^'^Ak) A2,sAk) \ , . 

^^P-^'^^- 1^ A2,«,.,(fc + Q) Ai,,,,,(fc + Q) )■ ^'"^ 

Although the off-diagonal matrix element of Aband(^) 

is finite in general, the low-energy properties below Tc 
are hardly affected by the off-diagonal components when 
|A| <C \a\. Therefore, we simply drop the off-diagonal 
components and obtain the Hamiltonian eq. (6). The SC 
order parameter for each band is expressed by the diag- 
onal components as, A^(fc) = (Aband(^))77- 
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The normal and anomalous Green functions are ex- 
pressed in the band basis as, 

Gj{k+,iUn) = {iu!n + ej{-k-))/A^{k,iLOn), (H) 

F^(fc,ia;„) = -A^(fc)/A^(fc,ia;„), (12) 

with 

Aj{k, iun) = (iw„ - e^{k+)){iu}n + ej{-k-)) - |A^(A:)p, 

(13) 

where u;„ = (2n + l)nT is the Matsubara frequency and 
T is the temperature. 

We decompose the uniform spin susceptibility into the 
Pauli part and Van-Vleck part, 



_ P I V 



(14) 



The Pauli susceptibility xji,^ arises from the intra-band 
scattering while the inter-band scattering gives rise to 
the Van-Vleck susceptibility (WS) xjiv- the following 
we assiimc the staggered moments along the principal 
axis, namely hq || x,y or z. Following the Appendix A, 
the Pauli susceptibility and VVS are expressed as, 



-limg-_^olimn„^o ^ 



7,fc 



X [G^{k + q)G^{k) ±F^{k + q)F^{k)], (15) 



and 



.=EE^-^^K^^) ^^^-fj^^^;f^^ (16) 

^^5 k es{k)-e^{k) 

respectively. The sign in cq.(15) is + for fi ^ x,y and — 
for fi = z.Wc define A':^'g{k) = S'^^g{k,k)S^^ik,k) where 
Sl^g{k + q, k) is the spin operator in the band basis. 



with 



S^'{k + q,k) = U\k + q)S''U{k), (17) 

(18) 



SI' = 



r(M) 











The expression of Pauli susceptibility eq. (15) is equiv- 
alent to the spin susceptibility in a multi-band system. 
The sign + and — in eq. (15) correspond to the cen- 
trosymmetric superconductor with spin singlet pairing 
and that with spin triplet pairing for d _L _ff, respec- 
tively. Thus, the Pauli susceptibility of NCSC decreases 
in the afe-plane below Tc while that is constant for the 
magnetic field along the c-axis. 

It should be noted that the VVS has a temperature 
dependence above Tc which is similar to that of Pauli 
susceptibility, because the ASOC is much smaller than 
the Fermi energy (|a| <C Ef)-^^ Therefore, the VVS in 
eq. (14) should be included in the spin part of mag- 
netic susceptibility which is extracted by the K-x plot. In 
this sense, the VVS, arising from the band splitting due 
to the ASOC, is quite different from the better known 
VVS coming from the orbital degrees of freedom. Note 
that both VVS are not affected by the superconductivity 
when Tc <C lal. 



If the order parameter is spatially uniform, namely 
gk = 0, the Pauli susceptibility is described by the mo- 
mentum dependent Yosida function as. 



7 



dk^A^^^Jk^)Y{^^{kF),T)/v^{k^), (19) 



for ij, = x,y, and 



Xz 



7 



dkFA^l^{kF)/v^{kF), 



(20) 



where / dfcp is the integral on the Fermi surface, and 
u^(fcF) is the Fermi velocity of 7-th band. The Yosida 
function is defined as, 



y(A,T) = -J d£/'(Vff2 + A2), 



(21) 



where f'{E) = df /dE is the derivative of the Fermi dis- 
tribution fimction. Since F(A, 0) = and F(0,T) = 1, 
we obtain X/i,i(r = 0) = for jj = x,y and Xzz(T = 
0) = X^z + XL{T = Tc) = Xzz(T = Tc). Thus, the resid- 
ual spin susceptibility along a&-plane is given by the VVS 
alone, while for fields parallel to the c-axis both the Pauli 
and Van-Vleck susceptibility contribute. It should be no- 
ticed that the spin susceptibility at T = is independent 
of the pairing symmetry. In this sense the spin suscepti- 
bility is universal in the linear response regime when the 
system lacks the inversion symmetry. 

2.2 PM state 

We concentrate now on the imiform state (qn = 0) to 
investigate the residual spin susceptibility X/i/i at T = 
for fj, = x,y. The helical SC state with ^ ^ will be 
discussed later in §8. In the PM state we set hq = 
and assign the four bands as ei_2(fc) = £{k) ± a|^(fc)|, 
S3,4{k) = s{k + Q) ± a\g{k + Q)\ so that we can express 
the unitary matrix as. 



U{k) = 



U2{k) 
U2{k + Q) 



(22) 



where 

^'^^^ ^ 71 ( .9x(fc) + i.9y(fc) -5x(fc) - igyik) ) '^^^^ 

with gf_i{k) — gfj.ik)/\g{k)\. The matrix element of spin 
operator is obtained as, Sii{k,k) = — S'22(A:, fc) = ff^(fc) 
and (fc. A:) = -S^{k,k) = g^ik + Q) for n = x,y 
while S^^{k,k) = 0. We then find xL = Xyy = /5/2 at 
T = Tc where p is the DOS in the normal state. Since 
Xxx = Xxx + Xxx = p + 0{a^ /sp) in the normal state, the 
residual spin susceptibility at T = is obtained as, 

Xxx(T = 0) = xL = Xxx(T = Tc)/2 + 0(aV4)- (24) 

Thus, the spin susceptibility in the a6-plane at T = 
is half of the normal state value in the limit |a| <C 
£f- Qualitatively the same result has been obtained in 
Refs. 18,20,22-28. Fujimoto has shown that the VVS in- 
creases when the DOS has strong asymmetry and |a| is 
moderate.^^ However, the /3-band of CePtsSi which we 
will investigate later does not satisfy this condition.^'' 
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Since the spin susceptibility decreases below for the 
magnetic field along the a6-plane, the paramagnetic de- 
pairing effect of Hc2 should be observed in NCSC with 
Rashba-typc spin-orbit coupling. This is consistent with 
the recent observation of the paramagnetic depairing ef- 
fect in CeRhSi3'^'^° and CelrSig^'*'^'^^ under high pres- 
sure where the AFM order is suppressed. However, this is 
not the case in CePtsSi at ambient pressure (within the 
AFM phase) This observation leads us to 
study the influence of AFM order in the following part. 

2.3 AFM state with inversion symmetry 

In order to clarify the influence of AFM order, we first 
investigate the spin susceptibility in the SC state with 
inversion symmetry for hq 0. Owing to the inversion 
symmetry, the residual spin susceptibility depends on the 
pairing symmetry in the usual way. Here we discuss the 
spin singlet pairing state while the spin susceptibility is 
constant through Tc in the spin triplet pairing state. The 
spin susceptibility consists of the Pauli part and Van- 
Vleck part as in §2.1, and the Pauli part vanishes in the 
spin singlet pairing state at T = 0. 

As a result of the simple calculation, we obtain for the 
Pauli susceptibility above Tc, 



= Y^[6{e,{k)) + S{e,{k))]=p, 



(25) 



for ^ II Kq and 



E 



[5(ei(fc)) + ^(e2(fc))], (26) 



e+(fc)±^/e_(fc)2 + /;,| with 



for H _L /iQ. Here, ei^2{k) 

Hq = \hQ\ and e±(fc) = (e(fc) ± e{k + Q))/2. The Pauh 
part of spin susceptibility for the magnetic field perpen- 
dicular to the AFM moment decreases with growing hq , 
i-e-, xl^ihQ = 0) > xl^^ihQ ^ 0) ioi H ± hq. A Van- 
Vleck part is induced by the AFM order and leads to 
the residual spin susceptibility for H _L Hq at T — 0. In 
contrast, the Van- Vleck part and the residual spin sus- 
ceptibility vanish for the magnetic field parallel to the 
AFM moment. 

In Fig. 1 we show the numerical results for the spin 
susceptibility in the SC state at T = 0. For this numerical 
analysis we assume a tight-binding model approximating 
the so-called /3-band of CePtsSi,^'^ 

e{k) = 2<i(cos kx + cos ky ) + 4i2 cos k^ cos ky 

-1-2^3 (cos 2kx + cos 2fcy) -I- [2t4 + 4t5(cos k^ + cos fcy) 

+4i6(cos 2kx + cos 2fcy)] cos kz + 2tr cos 2kz — jic- (27) 

The chemical potential /Xc is determined so that the elec- 
tron density per site is n. The Fermi surface of the 
band, which has been obtained in the band structure cal- 
culations without taking AFM order into account, 
is reproduced by choosing the parameters as 

{t\,t2,ts,t4,,t^,tQ,t'j,n) = 

(1, -0.15, -0.5, -0.3, -0.1, -0.09, -0.2, 1.75), (28) 



and a — 0.3 defining ti as the unit energy. As shown 
in Fig. 1, the spin susceptibility for the magnetic field 
perpendicular to the AFM moment m = Xxx((3, 0)/iq is 
increased in the SC state by the AFM order while that 
in the normal state is little affected. 



a=0 



0.32 




Fig. 1. Spin susceptibility along the b- and c-axcs at T = 
against the staggered magnetic moment m along the a-axis. The 
ASOC is zero (a = 0) and the staggered field hq = HqX is as- 
sumed. The staggered magnetic moment is defined as m = | < 

s' '^3s''^i s'^i,a' > I s° that m = 1 is the full moment. The 
solid and dashed lines show the spin susceptibility in the normal 
state and in the spin singlet SC state, respectively. The spin sus- 
ceptibility along the a^axis is zero in the spin singlet SC state at 
T = 0. 



At this point we can discuss the role of the band struc- 
ture. According to eq. (26), the Pauli part of the spin 
susceptibility is small for the magnetic field H ± hq, 
if the quasiparticle dispersion is quasi-two dimensional 
and £_(A:) is small. Although the /3-band of CePtsSi has 
a three dimensional Fermi surface, the band dispersion 
is weak along the fc^-axis according to the result of band 
calculation.^'^ This means that the AFM order signifi- 
cantly would affect the SC state in this band of CePtsSi. 
The presence of quasi-two dimensional Fermi surface is 
also expected in CeRhSis.^^ 

There are cases where AFM order plays indeed an 
important role in a centrosymmetric material. For ex- 
ample, UPd2Al3 is a spin singlet superconductor with 
Tc = 2K which coexists with the AFM order. xhe 
AFM state has a high Neel temperature of Tn = 14K 
and a large staggered magnetic moment, m = 0.85/LtB- 
This moment is directed to the afe-plane of tetragonal 
lattice and Q = (0, 0, tt).'''^ This is the same spin struc- 
ture as CePtsSi. NMR measurements show the decrease 
of Knight shift below Tc with a large residual part.*""* The 
VVS arising from the AFM order may induce the large 
residual spin susceptibility, although the multi-orbital ef- 
fect is another possible origin. This is consistent with 
the large iJc2 which exceeds the standard paramagnetic 
limit.s^ 

2.4 AFM state without inversion symmetry 

The result in §2.3 implies that the AFM order en- 
hances the Van- Vleck part of spin susceptibility in the 
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non-centrosymmctric system for the magnetic field per- 
pendicular to the AFM moment. We have shown the 
results for the spin susceptibility along the a- and 6- 
axes by assuming the dispersion relation eqs. (27), (28), 
a = 0.3 and Hq \\ x to describe the electronic structure of 
CePtaSi below Tn (Fig. 4 of Ref. 27). For fields H ±hQ 
the normal state and SC state susceptibility merge for in- 
creasing staggered moment, suggesting a diminishing of 
the reduction of the spin susceptibility in the SC state. 
On the other hand, the behavior is opposite for H \\hQ. 
Thus, a remarkable 2-fold anisotropy is expected in the 
spin susceptibility below even if the anisotropy is weak 
in the normal state. The condition H ± Kq is generally 
favored because the magnetization energy is maximally 
gained for the field direction with largest spin suscep- 
tibility. However, the meta-stable state H \\ hq can be 
realized in the weak magnetic field which is smaller than 
the anisotropy energy of AFM moment. 

The role of AFM order is suppressed by increasing the 
ASOC. We have confirmed that Xyy in the SC state is 
decreased by increasing a when Kq \\ x. The AFM order 
plays a quantitatively important role when the ASOC is 
much smaller than the Fermi energy. 

If the AFM moment is parallel to the c-axis as in 
CeCoGes,^ the spin susceptibility along both a- and 6- 
axes is increased in the SC state by the AFM order, while 
that along the c-axis is not affected. 

3. Effective Model for CePtsSi, CeRhSig and 
CelrSia 

In preparation for the discussion of the non-linear re- 
sponse to the magnetic field we will introduce here an 
eflfective model for CePtaSi, CeRhSia and CelrSia. This 
is important as we will show that the non-linear spin 
susceptibility significantly depends on the symmetry of 
order parameter in contrast to the universal spin suscep- 
tibility in the linear response theory (see §2.1). In partic- 
ular, we point out the strong non-linearity in the pairing 
state with dominant p-wave character. 

We analyze the following effective model, 

H = H^ + Hso + Haf + Hz + Hi, (29) 
Hz = -^h-Sik), (30) 

k 

Hi = uY^ni^fUi^i + {V- J/4) ^ mrij 

i <i,j> 

+J ^ (5^, • Sj - 2SfS^) (31) 

<i,J> 
k^k' ,q,s 

HU +{V- ^)C(fc - fc')}c^_,_,-4^/£;,c_,, 

where rij^s is the electron number at the site i with spin 
s, rii = rii^i +ni^i, s = —s and C{k) = 2(cosA;x + cosfcy). 
The spin operator in the real space basis is defined as 
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Si = J2s s' ^ss'cj gCi^s' ■ The bracket < i,j > denotes the 
summation for the nearest neighbor sites in the a6-plane, 
namely j = i±a ov j = i±b with a and b the unit vectors 
along the a- and &-axis, respectively. 

The first three terms in eq. (29) have been defined ear- 
lier in eqs. (2-4). For the dispersion relation, we adopt 
the tight-binding model eq. (27) with using the parame- 
ter set eq. (28) and a = 0.3, reproducing the /3-band of 
CePtaSi.''®"^^ We choose the /3-band, because it has sub- 
stantial Ce 4/-electron character^^ and the largest DOS 
at the Fermi energy, namely 70% of the total DOS.^° Be- 
sides the sizable jump in specific heat, also the remark- 
able isotropy of i?c2 between the a6-plane and c-axis®* 
also indicates that the three-dimensional Fermi surface 
of the /3-band is mainly responsible for the superconduc- 
tivity in CePtaSi. In Appendix B we will investigate the 
other dispersion relation which favors the <ix2_y2-wave 
superconductivity. 

As for the; AFM order, we assume the staggered field 
pointing along the [100]-direction Kq = hqx with Q = 
(0, 0, tt) following the experimental results of CePtaSi.®'' 
For the magnitude we assume Hq <C W, the band width. 
This is consistent with the small observed magnetic mo- 
ment ~ 0.16/ub in CePtaSi.®'^ The AFM moment is ex- 
pected to be small also in CeRhSia and CelrSia since 
superconductivity occurs near the AFM quantum criti- 
cal point. 

The fourth term Hz is the Zeeman coupling term 
due to the applied magnetic field. We have defined h = 
^qHbH where g is the ^-factor of quasiparticles and /^b 
is the Bohr magneton. The paramagnetic depairing effect 
on the superconductivity is characterized the dimen- 
sionless coupfing constant h/Tc with h= \h\. 

The last term Hi describes the effective interaction 
leading to the SC instability and includes three coupling 
constants, U, V and J. We assume the on-site interaction 
U and the interaction between the nearest neighbor sites 
in the a6-plane V. The coupling constant J describes the 
part of interaction arising from the AFM order which is 
anisotropic. According to the random phase approxima- 
tion (RPA) for the Hubbard model,^^''^^ the SC order 
parameter is affected by the AFM order mainly through 
the anisotropy of effective interaction, which can be de- 
scribed by the J-term in eq. (31). 

In the following we examine two parameter sets, 

(A) U>0, V = -0.8U, (33) 

(B) U<0, V = 0. (34) 

The amplitude of U is chosen so that Tc = O.Ol at zero 

magnetic field. The ground state is dominantly (A) p- 
wave and (B) s-wave, respectively. Hereafter we simply 
call these states p-wave and s-wave state, respectively. 
The parameter set (A) is the most important for our pur- 
pose, because the p- wave symmetry is the most promising 
candidate for the pairing state in CePtaSi.^^' 27.29.85, 86 
Although the spin triplet superconductivity is handi- 
capped due to the lack of inversion symmetry in non- 
centrosymmetric systems according to the Anderson's 
theorem, *7 the depairing effect arising from the ASOC 
vanishes (or is at least smallest) in the p-wave state 
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with d{k) II g{k).^^ This condition is not satisfied in 
the reahstic model, however the depairing effect due to 
the ASOC is almost avoided in the p-wave state with 
d{k) = —pyX+p^yP"^ Another parameter set (B) is inves- 
tigated as a typical model for the dominantly spin singlet 
pairing state. We will investigate the dominantly (ix2_y2- 
wave state in Appendix B and obtain qualitatively the 
same results as the s-wave state. 

Before analyzing the effective model in cq. (29), we 
comment on the RPA theory applied to the Hubbard 
model for the /3-band of CePt3Si.27.62 xhis theory leads 
to two possible pairing states due to spin fluctuation me- 
diated interaction: the ,s + _P-wave and the p+Z) + /-wave 
state. The former is dominated by the p-wave component 
and can be viewed as an intra-plane pairing state, while 
the latter is described by the inter-plane pairing dom- 
inated by the d^z- and dyz-wave components. Here we 
focus on the s + P-wave state whose order parameter 
is reproduced by assuming the parameter set (A) and 
J = O.W (J = 0) for hq = 0.125 {hq = 0) in eq. (29). 
On the other hand, the p + D + /-wave state is not re- 
alized in eq. (29) because the inter-plane interaction is 
neglected. It is expected that the paramagnetic proper- 
ties in the inter-plane d-wave state are qualitatively the 
same as those in the intra-plane s- and rf-wave states. 
The other characteristic properties of the p + D + /-wave 
state will be discussed in §6 and §7. 

To solve the effective model eq. (29), we apply the 
mean field theory and obtain the mean field equations 
as, 



Ai,,,,-(fc) = -T^{C/ + (1/ - ^)C{k - k')} 

n,k' 



xFi^s,s{k' .l^n), (35) 

A,,,,,(^) = -Tj2c{k-k'){VFi,,,s{k',L0n) 



- 2^i>s>s(^'''^")}- 



(36) 



The normal and anomalous Green functions 
Gi,s,s'{k' ,LOn), Fi^s,s'{k' ,u>n) in the spin basis are 
obtained by the Dyson-Gorkov equation. 



Gn{k+,U!ny 



Aspin(^) 



\^ A^„(fc) -GT(-fc_,-a;„)-i ) 

( G{k+,uJn) F{k,u;n) \_, ,„„^ 

F^KoJr.) -G^(-fc-,-a;„) ) 

where X{k) {X = G, F, Agpin) is the 4x4 matrix, 

X(k)=( ^^^^Ak)^ X2,sAk)^ \ (3g^ 
^ ' \ X2,s,s'{k + Q) Xi,,,,.(/e + Q) J ^ ' 

The normal Green function in the normal state, 
G'N(fc,w„) is obtained as G'N(fc, w„) = (iw„l — H{k))~^ 
by using eq. (8) with TiQ = hqx and e{k) = e{k)a^^^ + 
agik) ■ a — h ■ a. 

We here discuss the symmetry of the SC state on the 
basis of the following parameterization of order parame- 




-dx(fc) + irfy(fc) 

-$(fc) + 4(fc) 



<^{k) + d,ik) 

d^{k) + idy(fc) 



(39) 



where we use the even parity scalar function $(fc) and 
the odd parity vector d{k). Although a second component 
A2, s,s' {k) appears in the AFM state, the basic properties 
and symmetries are little affected by A2,s,s'(fc). 



i^-\va\x! si al 


E\'('ii i)aril\' i)arl 


Odd i>arii\' ]>arl 


PM at /i = 


k{S + Cx + Cy) 


(-Sy,Sx,0) 


AFM at ft = 


k{S + r)Cx -1- Cy) 


(-Sy,/3Sx,0) 


AFM at h = hy 


k{S + ?7Cx + Cy) 


(-Sy,/3Sx,-i7Sy) 



Table I. Symmetry of order parameter in the dominantly p-wave 
state. We use the abbreviations Cx,y = cos fex,y and Sx,y = 
sinfcx,y. We assume the PM state at zero magnetic field, the 
AFM state at zero magnetic field and the AFM state under the 
magnetic field h \\ y from the top to the bottom. The parameters 
/?, 7, K, 5 and 77 are real. 



In Table I we summarize the order parameters in the 
p-wave state. The admixture with the even-parity com- 
ponent due to the ASOC is expressed by the parameter 
K which is in the order of a/ep- We obtain k ~ 0.15 for 
a = 0.3. The even-parity part ^(k) is dominated by the 
extended s-wave component and (5 ~ 0.2 because the con- 
ventional s-wave component is suppressed by the strong 
on-site repulsion U. The dz-component of the odd parity 
vector d{k) is induced by the magnetic field 7 oc /i/ep to 
gain the Zeeman energy. The SC state is mainly affected 
by the parameter f3 which is unity in the absence of AFM 
order and magnetic field. The magnetic field along the 
[010]-axis ([100]-axis) decreases (increases) (3. The influ- 
ence of the AFM order depends on the value of J. We find 
P ~ 0.3 (P ~ 3.7) for J = 0.3y (J = 0) at hq = 0.125 
and h = 0. The deviation from /3 = 1 can be viewed 
as the mixing between d{k) = (— sin fcy , sin fcx, 0) and 
another p-wave state d(k) = (sin fcy, sin /cx, 0). Although 
these pairing states belong to different irreducible repre- 
sentations of the symmetry, they are mixed due to 
the presence of the symmetry reducing AFM moment or 
magnetic field. 

In general, the dominantly s-wave state is admixed 
to the p-wave state due to the ASOC and belongs to 
the same irreducible representation as the dominantly 
p-wave state realized for the parameter set (A). How- 
ever, only the conventional s-wave component ^{k) = 1 
appears and d{k) = for the parameter set (B) be- 
cause only the on-site interaction is taken into account 
{V = J = 0). Generally speaking, the admixture of spin 
singlet and triplet order parameters plays no important 
role when the ASOC is much smaller than the Fermi en- 
ergy, jaj < £f. 

The "helicity" ^-vector is perpendicular to the mag- 
netic field as will be discussed in §8 in details. The ampli- 
tude of qn below should be determined to maximize 
the condensation energy. However, here we determine qu 
at T = Tc{h) and neglect the temperature dependence 
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below Tc for simplicity. The transition temperature Tc{h) 
is determined by linearizing the mean field equation as, 



]{U+{V-^-)C{k-k')} 



n.k' 



(40) 



n.k' 



X{V(t)i^s,s{k' ,UJn) - —4>i,s,s{k',LLln)}, (41) 

where (f)i^s.s'{k' ,uJn) is obtained by linearizing the 
anomalous Green function Fi^s.s'{k' ,ujn) with respect to 
Aspm(fc')- We optimize the eigenvalue \{q) with respect 
to the order parameter Aspin(fc) and the hehcity g = gn- 
The transition temperature Tc{h) is determined by the 
criterion X{qH) = 1. 

We have estimated the condensation energy below Tc 
and found that the magnitude of (fn increases as decreas- 
ing the temperature. However, we have confirmed that 
the temperature dependence of qh can be ignored for 
the magnetic properties discussed in the following part. 

4. S-wave State 

For the discussion of non-linear response to the mag- 
netic field in NCSC we first discuss the simplest case, 
namely the s-wave state without AFM order. We ad- 
dress the enhancement of the critical magnetic field 
hc2 = 53/^6-^02 due to the ASOC, assuming the param- 
eter set (B) U <0 and V = 0. Figure 2 shows the phase 
diagram, temperature T/Tc versus magnetic field h/Tc 
along the [100]- or [010]-direction. The critical field hc2 
for both the uniform state (qh = 0) and the helical state 
((fn 7^ 0) are depicted, whereby also the behavior in the 
absence of ASOC {a = 0) is included for a comparison. 
As we focus here on the paramagnetic limiting effect, 
we neglect the orbital depairing for simplicity. Note that 
hc2 in the helical s-wave state has been investigated in 
Ref. 35 for an isotropic Fermi surface. 

The data in Fig. 2 demonstrate that the hc2 is signifi- 
cantly enhanced by the ASOC. This is partly due to the 
residual spin susceptibility in the SC state induced by 
the ASOC. Neglecting the magnetic field dependence of 
the spin susceptibility, we obtain a simple estimation for 
the critical magnetic field. 



hc2 = 



2Ec 



(42) 



where is the condensation energy and and x'^ 
are the spin susceptibility in the SC and normal state, 
respectively. According to eq. (42), hc2 increases by a 
factor of \/2 because of the residual spin susceptibility 
= ^X^ at T = 0. In fact, hc2 is enhanced even more 
due to the magnetic field dependence of spin suscepti- 
bility. A further enhancement of /ic2 is caused by the 
formation of a helical SC state, which exceeds the en- 
hancement in centrosymmetric superconductor owing to 
the presence of an FFLO state.®* 



a=0.3 uniform 
a=0.3 helical 
a=0 uniform 
a=0 FFLO 




Fig. 2. (Color online) The ff-T phase diagram in the s-wave state 
for the magnetic field along the [100]- or [010]-axis. The dia- 
monds (circles) show the reduced critical magnetic field hc2/Tc = 
^g^BHc2/Tc in the helical (uniform) SC state against the re- 
duced temperature T/Tc- We assume C/ < 0, V = 0, J = 0, 
a = 0.3 and /iq = 0. The phase diagrams in the absence of 
ASOC (o = 0) are shown for a comparison. The dashed and 
solid lines show the hc2/Tc in the uniform (gh = 0) and FFLO 
(& 7^ 0) states, respectively. 



o 4 



♦ -♦PM,h//[010] 
•-•AFM, h//[010] 

^ A AFM, h//[100] 




Fig. 3. (Color online) The H-T phase diagram in the helical s- 
wave state with AFM order. We assume /iq = 0.125x. The other 
parameters are the same as Fig. 2. The circles and triangles show 
the hc2 for the magnetic field along the [010] and [100]-axis, 
respectively. The hc2 in the PM state is shown for a comparison 
(diamonds). 



We here investigate the influence of AFIVI order on the 
s-wave SC state. Figure 3 shows that the hc2 in the s- 
wave state is increased by the AFIVI order, however the 
enhancement is very small. According to the simple es- 
timation eq. (42) and the universal spin susceptibility 
in the linear response theory (see §2.4), the hc2 along 
the [010]-axis ([100]-axis) is enhanced (suppressed) by 
the AFIVI order through the increase (decrease) of x^- 
However, the enhancement (suppression) is much smaller 
than expected in this simple estimation. This is mainly 
because of the formation of helical SC phase which in- 
duces the non-linear spin susceptibility at high fields. In 
§5.3 we will show that the infiuence of AFIVI order is 
much more significant in the dominantly p-wave state. 
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5. The p-wave State 

We here investigate the dominantly p-wave state which 
is the most promising candidate for the pairing state in 
CePtaSi. In this section we assume the parameter set (A) 
U >0,V = -0.8U. 

5.1 PM state 




T/T 

c 

Fig. 4. (Color online) The H-T phase diagram in the helical p- 
wave state (circles). We assume U > 0, V = — O.Sf/, J = 0, a = 
0.3 and hq = 0. The hc2 in the the s-wave state (diamonds) and 
in the d^2_y2-wave state (squares) are shown for a comparison. 
The parameter set for the d^2_y2-wave state is shown in the 
Appendix B. 



To illuminate the difference with the dominantly spin 
singlet pairing state we again turn to the PM state. We 
find that paramagnetic depairing effect is naturally less 
effective in suppressing the onset of superconductivity. 
Figure 4 shows hc2 for the p-wave state which is much 
higher than for the case of s-wave as well as d-wave pair- 
ings. This is rather surprising because the /ic2 is inde- 
pendent of the pairing symmetry considering only the 
simple estimation in eq. (42). Actually, hc2 of the p-wave 
state is enhanced by the modification of SC order param- 
eters due to the mixing with d{k) = (sin /cy, sin fcx, 0) in 
addition to the formation of helical SC state. Since the 
p-wave superconductivity has a multi-component order 
parameter with respect to the spin, the order parameter 
can be modified to optimally cope with the competition 
between the Zeeman coupling energy and ASOC. This is 
not the case in the dominantly spin singlet pairing state. 
This is the main reason why the paramagnetic depairing 
effect in NCSC depends on the symmetry of the leading 
order parameter. We see that the hc2 curves in Fig. 4 
merge in the low magnetic field region where the linear 
response theory is justified. 

In order to shed light on the mechanisms stabilizing 
the p-wave superconductivity at high magnetic fields, 
i.e., (i) the formation of helical SC state, and (ii) the 
modification of SC order parameters, we compare hc2 
with the one for the uniform state with (fn = (tri- 
angles in Fig. 5) and the one for the SC state with 
d{k) = (— sinA:y,sinA:x,0), <^{k) = and gk = (dia- 
monds in Fig. 5). Both (i) and (ii) are neglected in the 
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T/T 

c 

Fig. 5. (Color online) The circles (triangles) show the H-T phase 
diagram in the helical (uniform) p-wave state. We show the phase 
diagram in the SC state with d{k) = (— sin fcy , sin fcx, 0), ^{k) = 
and (jH = for a comparison. The parameters are the same as 
in Fig. 4. 



latter (diamonds) while (i) is neglected in the former (tri- 
angles). The comparison between the triangles and dia- 
monds shows the enhancement of hc2 by optimizing the 
SC order parameter. Actually, the d^- (dy-)component of 
d-vector decreases in the magnetic field along the [100]- 
([010]-)axis to avoid the paramagnetic depairing effect. 
The hc2 is furthermore enhanced below T = O.GTc by 
forming the helical SC state (circles). Thus, the p-wave 
superconductivity can be stabilized in the magnetic field 
which is much higher than the standard paramagnetic 
limit owing to the combination of mechanisms (i) and 

(ii)- 



1 1 1 1 1 1 1 1 1 1 


. ■ • ^ 




^ 


I.I.I 


— h=O.IT^ 

— - h=lT^ ' 

. -h=I.5T^ " 

— h=2T 

C 

1 , 1 



\ 0.2 0.4 0.6 0.8 I 

T/T 



c 

Fig. 6. (Color online) The spin susceptibility along the [100]- and 
[010]-directions in the helical p-wave SC state without AFM or- 
der. The magnetic field is chosen as /i = O.lTc, h = T^, h = 1.5Tc 
and h = 2Tc from the bottom to the top. The other parameters 
are the same as in Fig. 4. 



A large critical magnetic field hc2 generally indicates 
that a SC state with a large spin susceptibility is stabi- 
lized at high magnetic fields. The general spin suscepti- 
bility defined by Xa = Xb = My^^y/h is obtained from the 
calculation of the uniform magnetization, 

M^ = Y,TvS^'G{k), (43) 

k 
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where is the spin operator defined in eq. (18). The 
corresponding spin susceptibihty for the p-wave state is 
shown in Fig. 6. The spin susceptibihty ed h — O.lTc 
drops to half of its normal state value at T = 0. This is 
consistent with the linear response theory in §2.2. For a 
magnetic field comparable to or higher than the standard 
paramagnetic limit, the order parameter of p-wave state 
is modified in order to avoid the paramagnetic depairing 
effect. Therefore, the spin susceptibility at /i = 2Tc is al- 
most constant through T^., although the critical tempera- 
ture remains high, (Tc(ft. = 2Tc) = 0.61Tc{h = 0)). These 
results should be contrasted to the dx2_y2-wave case dis- 
cussed in Appendix B. The Tc of dx2_y2-wave state is 
reduced more strongly (Tc(/i = 2Tc) = 0A2Tc{h = 0)), 
but the decrease of spin susceptibility below Tc is larger 
than that in the p-wave state (Fig. B.l). 

5.2 Role of anisotropic Fermi surface 

A further important role in this context is played 
by the shape of the Fermi surface. The band structure 
of the /3-band is complicated but has one eye-catching 
property: the cross sections of the Fermi surface in the 
range > 27r/3 are quadrilateral with the corners 
along the [110]- and [l-10]-directions.^^^^° We show the 
schematic figures for the Fermi surface in Fig. 7 where 
the anisotropy is stressed for simplicity. The anisotropy 
of the Fermi surface affects hc2 in two ways, by facili- 
tating (i) the formation of the helical SC phase and (ii) 
the modification of the SC order parameters, as we will 
discuss now. 



(a) H=0 



(b) H//[010] (Uniform) 
I 





(c) H//[010] (Helical) (d) H//[110] 





Fig. 7. (Color online) The schematic figure for the Fermi surface 
and SC gap in the p-wave state. We show the cross section for a 
fixed fcz . Two solid lines show the Fermi surfaces which are split 
by the ASOC. The dashed lines show the magnitude of SC gap on 
each Fermi surface, (a) The uniform BCS state at _ff = 0. The 
direction of d-vector is shown by the arrows, (b) The uniform 
BCS state for _ff//[010]. The parts of Fermi surface "F and "IF 
are shown. The SC gap on the part "IF is suppressed, (c) The 
helical SC state for H//[010]. The SC gap on the part "II" of 
large Fermi surface is increased, (d) The helical SC state for 
H//[1W]. 



First, (i) the helical SC phase is stable for the 
anisotropic Fermi surface not only in the p-wave state 
but also in the s- and d-wave states. This is simply be- 
cause a set of quasi-particles with k — ±k + qji/2 can 
have low energy on a large part of the first Brillouin 
zone (nesting feature of the Fermi surface). As shown in 
Fig. 7(b), the large (small) Fermi surface moves to the 
right (left) in the magnetic field along the [010]-axis. Un- 
der this condition, uniform Cooper pairing on the Fermi 
surface part parallel to [010]-axis (part "11" in Fig. 7(b)) 
is destabilized, while it is little affected on the part "I". 
However, the depairing effect arising from the large Fermi 
surface is essentially avoided in the helical SC phase hav- 
ing qji ^ 2h/vpx (Fig. 7(c)) because of the nesting of 
Fermi surface along the [100]-direction. This leads to the 
strong enhancement of hc2- This is not the case in the 
isotropic system where the Fermi surface is not nested. 

Second, the anisotropic Fermi surface enhances (ii) the 
mixing of order parameters and increases in this way hc2 
in the p-wave state. Because of the structure of vector 
g{k) (X {—Vy{k),Vx{k),0), the SC gap on the Fermi sur- 
face perpendicular to the [010]-axis (Fermi surface "I" 
in Fig. 7(b)) is mainly induced by the dx-component 
of spin triplet order parameter while the dy-component 
is the main source of the SC gap on the other part 
(Fermi surface "II" in Fig. 7(b)). Since the d^- and dy- 
components induce the Cooper pairing on different parts 
of the Fermi surface, the coupling is weak between these 
two order parameters. Hence, the splitting of energy be- 
tween d{k) — (— sinfcy, sinfcx, 0) (most stable state) and 
d{k) = (sinfcy, sin fcx, 0) (second most stable state) due 
to the ASOC is small, and they can be easily mixed by 
the applied magnetic field. 



Isotropic FS 



Tetragonal FS 



d(l<) = (kx,ky,0). 



(ky.kx.O) . 
(-kx,ky,0) 



(-ky,kx,0) . 



(sinkx,sinky,0) 

(-sinkx,sinky,0) 
(sinky,sinkx,0) 

(-sinky,sinkx,0) 



Fig. 8. (Color online) The schematic figure for the energy lev- 
els in the dominantly p-wave state. The isotropic and tetrag- 
onal symmetries are assumed in the left and right figures, re- 
spectively. The 2-fold degeneracy in the isotropic system be- 
tween d{k) = (fcyjfcxjO) and d{k) = (— fcx,A:y,0) is lifted to 
d{k) = (sin fcy , sin fcx , 0) and d(k) = (— sin fcx, sin fey, 0) in the 
tetragonal system. In case of the /3-band of CePtsSi, d{k) = 
(sin fey , sin fcx , 0) has lower energy. 



In general, the tetragonal anisotropy of the Fermi sur- 
face reduces the splitting between the most stable and 
the second most stable pairing states. For an isotropic 
Fermi surface, the second most stable pairing state 
has 2-fold degeneracy; d{k) = (fcy,A:x,0) is degener- 
ate with d{k) = (— fcx,fcy,0). However, this degeneracy 
is lifted by the tetragonal anisotropy as shown in the 
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schematic figure (Fig. 8). This lift of degeneracy de- 
creases the difference of condensation energy between 
d(k) — (— sinfcy, sin fcx, 0) and d{k) — (sin /cy, sin /cx, 0) 
(or d{k) — (— sin fcx, sin A:y , 0)) in the system with tetrag- 
onal symmetry (see Fig. 8). 




T/T 

c 

Fig. 9. (Color online) The critical magnetic field hc2 along the 
[110]-axis in the helical p-wave (open circles) and s-wave (open 
squares) states. Those along the [100]-axis are shown by the 
closed symbols for a comparison. 



Furthermore, a strong anisotropy of the Fermi surface 
induces a pronounced 4-fold anisotropy in the paramag- 
netic properties. Figure 9 shows that the /ic2 along the 
[110]-direction is much smaller than that along the [100]- 
direction in case of the p-wave state. This is mainly be- 
cause the state d{k) = (— sin/cy, sin/cx, 0) is admixed by 
the magnetic field along the [110]-direction, with d{k) — 
(— sin fcx, sinfcy, 0), but not with d{k) = (sin fcy, sin fcx, 0). 
The latter is less stable than the former in case of 
the /3-band of CePtaSi. On the other hand, the 4-fold 
anisotropy is weak for the s-wave state as shown by the 
squares in Fig. 9. This indicates that the anisotropic 
Fermi surface enhances the hc2 in the p-wave state mainly 
through the mixing of SC order parameters. 

Finally, we comment on orbital depairing which we 
have neglected so far. The orbital depairing effect is 
reduced by the mixing of order parameters in the p- 
wave state. For example, the parameter /3 in d{k) = 
(— sinfcy, /Ssinfcx, — i7sinfcy) is decreased by the mag- 
netic field along the [010]-direction, and reduces the or- 
bital depairing effect, because the coherence length along 
the [100]-direction shrinks. Thus, the hc2 in the p-wave 
state is enhanced by modifying the order parameter 
through the suppression of the orbital depairing effect 
as well as the paramagnetic depairing effect. 

5.3 AFM state 

In the discussion of the influence of AFM order on the 
p-wave SC state we focus on staggered moments along 
the [100]-axis with the magnetic field parallel to the [010]- 
axis, since the AFM moment favors to be perpendicular 
to the field. The situation of the magnetic field parallel 
to the moment is described in §7. 

The influence of the AFM order on the p-wave state 
significantly depends on the anisotropic spin-spin inter- 



action, the J-term in eq. (31). The critical field /ic2 de- 
picted in Fig. 10 with hq = 0.125 in the AFM ordered 
phase shows a clear trend. While the AFM leads to a 
reduction of hc2 in the absence of the anisotropic spin 
interaction (J = 0), a strong enhancement is obtained 
for J = 0.3. 



6 



1 • ' 1 


1 1 


•— •PM 

■— ■ AFM, J=0.3V 
AFM, J=0 





'0 0.2 0.4 0.6 0.8 



T/T 

c 

Fig. 10. (Color online) The H-T phase diagram in the p-wave 
state for the magnetic field along the [010]-axis in the presence 
of AFM moment along the [100]-axis. The squares and triangles 
show the hc2 for J = 0.3V and J = 0, respectively. We fix 
hQ = 0.125 and choose the other parameters as in Fig. 4. The 
hc2 in the PM state is shown for a comparison (circles). 



We understand these results by analyzing the param- 
eter P of d{k) = (— sinfcy, /3sinfcx, 0) at zero magnetic 
field. 

For (3 < 1, the superconductivity is dominant on Fermi 
surface region "I" in Fig. 7(b), while the magnetic field 
along the [010]-axis suppresses Cooper pairing on the 
Fermi surface "11" . For this reason, this SC state is robust 
against the magnetic field along the [010]-axis. The mag- 
netic field reduces /3 even more enhancing the anisotropy 
of the SC gap. The enhancement of hc2 due to the AFM 
order is much more significant than expected in the sim- 
ple estimation eq. (42). In fact, the suppression of para- 
magnetic depairing effect in case of /3 < 1 can be viewed 
as a result of the strong non-linear response to the mag- 
netic field. The small energy scale (3Tc appears in this 
case and induces the strong non-linearity. This is the 
reason why the influence of AFM order is much more 
important in the p-wave state than in the s-wave state. 
If we assume J — 0, the parameter (3 is more than unity, 
which is incompatible with our RPA analysis for the Hub- 
bard model. On the other hand, we obtain /3 ~ 0.3 for 
J = 0.3V and hq = 0.125, giving the result consistent 
with the RPA theory. 

The strong non-linear response to the magnetic field 
clearly appears in the magnetic field dependence of spin 
susceptibility. We show the spin susceptibility Xb for J = 
0.3V {13 < 1) and J = (/3 > 1) in Figs. 11 (a) and 
(b), respectively. For low magnetic fields (h = O.lTc) 
Xb is enhanced by AFM order in both cases consistent 
with the linear response theory (§2.4). We find that Xh 
is furthermore enhanced for the moderate magnetic field 
h ^ Tc with (3 < 1 (Fig. 11(a)) although the critical 
temperature rc(^) is little decreased. According to these 
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Fig. 11. (Color online) The spin susceptibility along the [010]- 
direction in the helical p-wave state with AFM order. We assume 
J = 0.3V leading to /3 < 1 in (a) and J = leading to /3 > 1 in 
(b), respectively. The magnetic field is shown in the figures. The 
other parameters are the same as in Fig. 10. 



theoretical results, the NMR Knight shift measurement^^ 
by Yogi et al. and the /iSR measurement^^ by Higemoto 
et al. were carried out in the non-linear response regime. 

In contrast to /3 < 1, the moderate magnetic field h = 
Tc little affects the spin susceptibility (Fig. 11(b)). The 
non-linearity of spin susceptibility appears only in the 
high field region close to the critical magnetic field. This 
is a characteristic property of the SC state with strong 
paramagnetic depairing effect such as the spin singlet 
pairing state in centrosymmetric system. Qualitatively 
the same magnetic field dependence is obtained in the 
dominantly dx^-ya-wave state (see Fig. B.l in Appendix 
B). 

6. Pairing symmetry in CePtaSi 

Measurements of i?c2 and the Knight shift are consis- 
tent with p-wave superconductivity in CePtaSi at ambi- 
ent pressure. The temperature dependence of iJc2^'^'^'^* 
implies the absence of paramagnetic depairing. NMR and 
^SK Knight shift data show no decrease below Xi;^66,67 
though Tc remains rather high at applied magnetic fields. 
These findings could be understood based on the p-wave 
state with AFM order for which the theoretical results 
have been shown in Figs. 10 and 11(a). 

We here note that the other possible mechanisms 
for the high critical field Hc2 are unlikely relevant in 



CePtsSi. For example, a small g-factor has been sug- 
gested for CeCoIns (g ^ 0.63).^^ It is expected that the 
(7-factor of CeCoIus is significantly renormalized by the 
strong AFM correlation in the a&-plane.^" However, this 
is not the case in CePtaSi where the spin correlation 
in the a&-plane are dominantly ferromag netic.2'^'5'^ The 
strong coupling effect which has been ignored in this 
paper is another possible cause of high Hc2- But, the 
jump of the specific heat at T = Tc does not indicate 
strong coupling effects in CePtaSi,^^^'^'^^ in contrast to 
CelrSia.^i 

The p-wave state is consistent with the coherence peak 
in NMR i/r^T2i'29,92 ^^d the line node behaviors in 
various quantities.^^''^"''^^'^^"^^ Moreover the microscopic 
theory within an RPA theory suggests an in-plane p-wave 
state induced by the /3-band of CePtaSi.^^ 

For _L /iQ the experimental magnetic properties 
of CePtaSi at ambient pressure are consistent with the 
p-wave state with /3 < 1. This indicates the strong 
anisotropy of the effective spin interaction, which is de- 
scribed by the J-term in eq. (31) and is compatible with 
the RPA analysis. It does however not agree with the 
naive second order perturbation theory which leads to 
the p-wave state with /3 > 1.®^ This is because the role 
of spin fluctuation is underestimated within the pertur- 
bation theory.^'^ Based on this fact wc may state that 
there is some evidence for spin-fluctuation-mediated su- 
perconductivity in CePtaSi. 

When the magnetic field is parallel to the AFM mo- 
ment H II /iQ, the paramagnetic depairing effect is en- 
hanced (suppressed) in the p-wave state with /3 < 1 
(/? > 1). We have confirmed that the hc2 for H \\ hq with 
/3 > 1 is qualitatively the same as the hc2 for H J- hq 
with (3 < 1 (squares in Fig. 10). If the sample had a do- 
main structure with respect to the direction of AFM mo- 
ment, the SC state with the maximal Tc would mark the 
SC transition. Under such circumstances p-wave states 
with both (3 > 1 and /3 < 1 could "avoid" the paramag- 
netic depairing effect and would be consistent with the 
experimental results in CePtaSi.^'^''^'^^'^^'^^ 

We here comment on the inter-plane d-wave state 
which we found as another possible pairing state on the 
basis of the RPA theory.^^ Although the 2-fold degen- 
eracy exists in this state (dxz- and dyz-wave), the order 
parameter has no internal degree of freedom with respect 
to the spin. Therefore, the paramagnetic depairing effect 
cannot be avoided by modifying the order parameter in 
contrast to the p-wave state. Hence, the magnetic prop- 
erties are qualitatively the same as those in the s-wave 
state which seem to be incompatible with the experimen- 
tal results in CePtaSi. The inter-plane d-wave state is in- 
compatible with the coherence peak in the NMR 1/TiT 
too.^^ 

7. Proposals for test experiments 

Here we discuss several experiments which could help 
to establish the pairing symmetry for CePtaSi as well as 
CeRhSia and CelrSia. 

The influence of antiferromagnetism on the magnetic 
properties can be tested by using the fact that AFM 
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order can be suppressed by pressure in these materi- 
als.^"^' It follows from our results in §4, §5.1 and 
Appendix B, that in the purely SC phase the paramag- 
netic depairing should limit the upper critical field for 
H \\ ab and the spin susceptibility should decrease be- 
low Tc in the low-magnetic field regime. Actually recent 
measurements of i?c2 along a6-plane in CeRhSia^'^" and 
CelrSia'''^*'^® imply a clear paramagnetic depairing ef- 
fect in the purely SC region, consistent with the theo- 
retical view. Notably paramagnetic depairing seems less 
effective in the AFM state of CelrSia.^^ This is compat- 
ible with p-wave pairing. No studies of this kind have 
been performed so far for CePtaSi. 

A further aspect is the 2-fold in-plane anisotropy in 
the AFM state. Since the [100]- and [010]-axes are not 
equivalent in the AFM state, a 2-fold anisotropy appears 
in the a&-plane. Although the AFM moment perpendicu- 
lar to the uniform magnetic field is generally favored, the 
situation H \\ can nevertheless be realized for mag- 
netic fields low enough to leave the orientation of the 
AFM moment unchanged. 

We summarize the 2-fold anisotropy expected for each 
pairing state in Fig. 12 taking also the orbital depairing 
effect into account. We assumed here that the Hc2 de- 
termined by the orbital depairing is much higher than 
the standard paramagnetic limit field in CePtsSi,^'^'^'^* 
CeRhSig'^'^ and CelrSig.^'^^ The upper critical field due 
to orbital depairing is naturally enhanced by the heavy 
mass of quasi-particlcs in these heavy Fermion com- 
pounds. Under such conditions paramagnetic depairing 
can play a role in the high-field regime. 

Fig. 12(a) shows the H-T phase diagram in the p-wave 
state. For H \\ Hq the paramagnetic depairing effect is 
enhanced (suppressed) with /3 < 1 (/? > 1). Note that 
the opposite occurs for H J- hq (see Fig. 10). There- 
fore, a significant 2-fold anisotropy of Hc2 could appear 
at high magnetic fields for either /3 < 1 or /3 > 1, pro- 
vided the AFM moment remains pinned. Qualitatively 
the same anisotropy would occur at low magnetic fields, 
because the orbital depairing effect is anisotropic owing 
to the in-plane anisotropy of coherence length, namely 
the difference of and ^b- On the basis of the RPA the- 
ory for CePtaSi^^ we have estimated the anisotropy as 
dH,2/dT\T=TAH II [100]) : dH,^/ dT\T=TAH \\ [010]) = 

: = 0.672 : 1 at /iQ = 0.125. Thus, the H-T phase 
diagram is highly anisotropic in both high and low mag- 
netic field region as shown in Fig. 12(a). 

The strong 2-fold anisotropy in the a6-plane appears 
also in the inter-plane d-wave state due to the anisotropy 
of the coherence length. The 2-fold degeneracy between 
the dxz" and dy^-wave states is lifted by the AFM or- 
der. The staggered moment along the [100]-axis favors 
the dxz-wave state and yields a coherence length which 
is longer along the [100]-axis than along the [010]-axis. 
For this reason iJc2 close to T = Tc is smaller for the 
magnetic field along the [010]-axis. This anisotropy is 
suppressed at high magnetic fields because the paramag- 
netic depairing effect is nearly isotropic as in the s-wave 
state (Fig. 3). These considerations lead to the schematic 
phase diagram in Fig. 12(b). 



(a) p-wave (P<1) 



(b) d„-wave 
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Fig. 12. (Color online) Schematic figure for the 2-fold anisotropy 
in the H-T phase diagram. We assume the AFM order along 
the [100]-axis. (a) The p-wave state with /3 < 1. The opposite 
anisotropy is expected for /3 > 1. (b) The inter-plane d-wave 
(lixz-wave) state, (c) Possible multiple phase transitions in the 
inter-plane d-wave state for H \\ [010]. (d) The intra-plane d-wave 
(i^x^— y2-wave or dxy-wave) and s-wave states. 



It should be noted that the in-plane anisotropy of Hc2 
in the inter-plane d-wave state does not vanish if the 
quantum critical point of the AFM order is approached. 
This is in contrast to the p-wave state where the in- 
plane anisotropy is suppressed by decreasing the AFM 
moment. In the vicinity of AFM quantum critical point, 
multiple phase transitions can occur for the inter-plane d- 
wave state as discussed in Ref. 27. These multiple phases 
in the H-T plane are shown in Fig. 12(c) for the magnetic 
field along the [010]-axis. Pure dxz- and dyz-wave states 
appear in the high-temperature region and in the high- 
magnetic field region, respectively. The chiral dxz ± idyz- 
wave state is stabilized at low temperatures and fields. If 
the multiple phase transitions were observed in the H-T 
plane or in the P-T plane, it would be a strong evidence 
for the inter-plane d-wave state. Although some indica- 
tions for a second SC transition have been reported in 
CePtaSi,^'*'^'*"^^ it remains unclear whether it represents 
an intrinsic property or is caused by the sample inhomo- 
geneity 

In contrast to the p-wave and inter-plane d-wave 
states, the 2-fold anisotropy of Hc2 is very weak in the 
intra-plane d-wave and s-wave states because the para- 
magnetic depairing effect as well as the orbital depairing 
effect are nearly isotropic. Therefore, we obtain a simple 
phase diagram in Fig. 12(d). 

Since the 2-fold anisotropy of Hc2 is quite different 
between the dominantly p-wave, inter-plane d-wave and 
intra-plane spin singlet pairing states, the future exper- 
iment in the AFM state could identify the pairing sym- 
metry in CePtsSi, CeRhSia and CelrSia. It should be no- 
ticed that this experiment can be performed in CePtaSi 
without applying the pressure. 
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8. Helical Superconductivity 

In this section we discuss the nature of the hehcal SC 
state which is a novel SC phase specific to NCSC. The SC 
phase with a finite total momentum of Cooper pairs is 
stabilized in the presence of Rashba-typc spin-orbit cou- 
pling under a magnetic field in the a6-plane."^^"^^ This 
state bears some similarity with the FFLO state^^ in 
centrosymmctric superconductors, but has also impor- 
tant differences. First, the helical SC phase is stabilized 
immediately above Hci which is much lower than Hc2 
in extremely type II superconductors. This is in contrast 
to the FFLO state which appears in a narrow region 
near Hc2 only. Second, the phase of SC order parameter 
is modulated as A(r) = Ae"?*^'" in the helical SC state 
(which is the same form as in the Fulde-Ferrel (FF) state) 
while the Larkin-Ovchinnikov (LO) state with the spa- 
tial modulation of the amplitude, A(r) = Acos^fn^ = 
A(e'&^"*-F e-'«H*')/2, is more stable than the FF state.^^ 
Because the two momenta qu and —qu are equivalent in 
the centrosymmctric system, the order parameter has a 
double q structure in the LO state. On the other hand, 
5h is not equivalent to —qu in the non-centrosymmetric 
system under a magnetic field. For this reason the helical 
SC phase is realized in the NCSC at least just below the 
critical temperature. At higher fields and low tempera- 
ture also a "stripe SC state" can be realized,^^' which 
is similar to the LO state. 

Experimental evidence for the FFLO state has been 
obtained for CeCoIns,^*^ more than forty years after 
the theoretical proposal.*^' This is partly because the 
FFLO state is suppressed by a weak disorder. The 
stripe SC state, which is resembles the FFLO state, can 
be suppressed by a weak disorder too. In contrast to these 
states the helical SC state is realized even in the disor- 
dered material, if the superconductivity is present. Al- 
though there is no experimental verification of the helical 
SC state in NCSC so far, the existence of the helical SC 
phase is a mandatory features from a theoretical point 
of view. 

Now we turn to the effect of finite qu on the param- 
agnetic properties. Although the influence of the helical 
superconductivity has been taken into account in §4 and 
§5, the following discussion will be important for a deeper 
understanding. 

One of the characteristic properties in the helical SC 
state is the presence of a finite spin magnetization. In the 
low magnetic field region this magnetization is expressed 
as M = Mo + x'H with finite Mq. For simplicity, we 
here consider the PM state and assume the SC order 
parameter without gap nodes. Then, the magnetization 
is obtained as, 

^0 = lY.9ik){Bi{k) ■ m - Mk) ■ qu) (44) 

k 

^lij dkFgikF)ivi{kF) ■ qu)/vi{kF) - (1 ^ 2)] (45) 
= D{z X qu), (46) 

with D oc a. We define B^(k) = d{e^{k)/ E^{k))/dk 
where 7 is a band index. As shown in eq. (46), the mag- 
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netization is oriented along the direction perpendicular 
to qu- 

The helical superconductivity also affects the differen- 
tial spin susceptibility = dM^/dH^ when the SC 
gap has a node. According to eqs. (11-13), the quasi- 
particles suffer a Doppler shift in the helical SC state 
and the single particle excitation energy is expressed as 

\Je^{k)'^ + |A^(fc)|2±e,(fc)-gH/2. Following eq. (15), the 
Pauli part of differential spin susceptibility is obtained 
as, 




X Yh (zi, (fcp ) , I A^ (Af ) I , T) 7^;^ (fcp ) , (47) 

for n = x,y where Yu{v, A, T) is the generalized Yosida 
function, 

Yu{v,A,T) = -^ J de[f'{Ve^ + A^ + v-qu/2) 

+/'(Ve2 + A2-,;.5H/2)].(48) 

Since Yii{v,A,0) = i^l -4A2/(i/- qn)^ for |A| < 

\v ■ (fH|/2, the Doppler shift boosts the differential spin 

susceptibility in the SC state with a gap node like in 
CePt3Si.63-65 

In fact, the uniform BCS state is favored at i? = and 

the helical SC state is induced by an infinitesimal mag- 
netic field owing to the linear coupling between the mag- 
netization and the helicity qa (eq. (46)) with qu _L H. 
Since the amplitude of qu is linear in the magnetic field 
\H\, the formation of helical SC state leads to a correction 
to the linear response theory in §2 of the uniform state. 
However, the correction is negligible when \a\ <C ep be- 
cause the amplitude of qu is small, \qu\ ~ {a/eF)h/vF in 
linear order of small parameter a/ep- 

The helicity can play a quantitatively more important 
role in the non-linear response regime, because the am- 
plitude of increases from Igkl ~ {a/sF)h/vF in the low 
field region to \qu\ ^ h/vp in the high field region with 
a rapid crossover around h ~ T^. For example. Fig. 13 
shows the magnetic field dependence of the helicity |^| 
in the p-wave state, with a sharp increase of the helicity 
above /i = Tc. As a result the critical field hc2 is signif- 
icantly enhanced at high fields as shown in Figs. 2 and 
5. 

The nature of the crossover from |^h| ~ {a/eF)h/vF 
to \qu\ ~ /i/wp becomes obvious viewing the momentum 
dependence of eigenvalues \{q) in eqs. (40) and (41). Fig- 
ures 14(a) and (b) show the numerical results in the PM 
and AFM states, respectively. In Fig. 14(a), A(^ pos- 
sesses a crossover from a single to a double peak struc- 
ture, yielding a rapid increase of the helicity. This result 
implies that the nature of the helical SC phase is differ- 
ent below and above the crossover magnetic field. Actu- 
ally, the "stripe SC state" can be stabilized above the 
crossover field. As shown in Fig. 14(b), the crossover 
from the single to the double peak structure is suppressed 
by the AFM order. The eigenvalue \{q) has a single peak 
even in the magnetic field much higher than the standard 
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Fig. 13. (Color online) The amplitude of the helicity \qYi\ just 
below Tc. The circles and squares show the results in the PM and 
AFM states, respectively. We assume the magnetic field along 
the [010]-direction which leads to gn along the [100]-direction. 
J = O.Sy is assumed in the AFM state. The other parameters 
in the PM and AFM states are the same as in Figs. 4 and 10, 
respectively. 



paramagnetic limit. This is simply because the AFM or- 
der suppresses the paramagnetic depairing effect in the 
p-wave state. 

We would like to point here that CePtaSi is a good 
candidate for an experimental observation of the heli- 
cal SC phase. Actually, the large critical field i?c2 leads 
to the helical SC phase with large (fn (l&l ~ hjvp) in 
a large part of the H-T phase diagram. It seems to be 
difficult to detect the helical SC phase with small helic- 
ity |9h| ~ (oL I £-p)h I v-p because the wave length is much 
longer than the coherence length. Thus the high field 
phase with IqhI ~ hjv-p is more promising for the ex- 
perimental observation. The high field phase is stable in 
the p-wave state above /i = Tc as shown in Figs. 4, 10 
and 13. However, this phase shrinks in the SC state with 
dominantly spin singlet pairing and/or the strong orbital 
depairing effect which leads to small 7?c2- 

9. Summary and Discussions 

We have investigated the paramagnetic properties in 
NCSC. The SC states with leading p-wave, d-wave or s- 
wave order parameter have been examined in view of the 
heavy Fermion superconductors, CePtaSi, CeRhSia and 
CelrSis. 

First, the linear response to the magnetic field has been 
investigated with the particular interest on the role of 
AFM order. The spin susceptibility is universal in the 
sense that it is independent of the pairing symmetry at 
T = 0, if the ASOC is much larger than the SC gap, 
and results from the band splitting due to the ASOC. 
The spin susceptibility below Tc is increased in the AFM 
state due to the folding of unit cell, if the magnetic field is 
applied perpendicular to the AFM moment. The result 
is opposite for the magnetic field parallel to the AFM 
moment. 

Second, we have shown that the non-linear response 
to the magnetic field depends on the symmetry of the 
leading SC order parameter. In particular, the spin sus- 
ceptibility and iJc2 for the p-wave state are significantly 
enhanced by non-linear effects through (i) the formation 
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Fig. 14. (Color online) The momentum dependence of the eigen- 
value A(g^ in the linearized mean field equation (eqs. (40) and 
(41)) with q = (i}x,0, 0). (a) the PM state and (b) the AFM 
state. The parameters are the same as in Fig. 13. 



of a helical SC state and (ii) the mixing of SC order pa- 
rameters. The anisotropy of Fermi surface can increase 
these non-linear effects, and strengthen the influence of 
AFM order. Taking these aspects into account, the exper- 
imental results^^^'^'^^^^^'^^ for CePtaSi at ambient pres- 
sure (within the AFM state) are consistent with the p- 
wave state admixed with a secondary s-wave component. 
This is the pairing state which has been proposed re- 
cently by Frigeri et al?'^ and identified by the microscopic 
RPA theory.^^'^^ It has been shown that this p-wave 
state is consistent with the line node behavior^^' ^^"^^ 
and also with the coherence peak in NMR l/TiT.2i'29,92 
Although the RPA theory has identified the inter-plane 
d-wave (dxz- and dyz-wave) state as further candidate, 
this state seems to be incompatible with the Knight shift, 
Hc2 and NMR 1/TiT measurements at ambient pressure. 
According to these comparisons between the theory and 
experiment, CePtsSi is rather likely the first identified 
spin triplet superconductor in Ce-based heavy fermion 
systems. 

We have proposed several experiments which can pro- 
vide further evidences for the pairing state in CePtaSi 
as well as in CeRhSia and CelrSia- The first proposal 
is the pressure dependence in various quantities. If the 
AFM order is a major cause of the unusual properties in 
CePtaSi, a pronounced pressure dependence is expected 
in NMR, specific heat, thermal transport, superfluid den- 
sity and so on. If CePtaSi has a leading p-wave order 
parameter, the following behaviors are expected above 
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the critical pressure P ~ 0.6GPa. (a) The Knight shift 
decreases below Tc for the magnetic field along the aft- 
plane and below the standard paramagnetic limit, (b) 
The paramagnetic depairing effect is enhanced for ^ || ab 
but not for H \\c. (c) The low-energy excitations due to 
the accidental line nodes are decreased. (d) The coher- 
ence peak in NMR l/T^T is enhanced by the isotropic 
SC gap.^^ These pressure dependences are not expected 
in the intra-plane d-wave (dx^-y^- and dxy-wave) and s- 
wave states. The pressure dependence (c) is expected also 
in the intcr-planc rf-wavc state and then the additional 
phase transition occurs in the P-T and H-T plane. 

Another proposal for a future experiment is the 2-fold 
anisotropy arising from the AFM order. The strong 2- 
fold anisotropy is expected in the dominantly p-wave 
state while the anisotropy is negligible in the s-wave and 
intra-plane rf-wavc states. In the inter-plane d-wave state 
the strong 2-fold anisotropy is expected near T = 
but the anisotropy is suppressed at high magnetic fields. 
The experimental observation of the 2-fold anisotropy in 
the AFM state could provide an important evidence for 
the pairing symmetry in CePtsSi, CeRhSis and CelrSis. 
Thus, the response to the AFM order can be a signature 
of the pairing symmetry in non-centrosymmetric super- 
conductors. 

Acknowledgments 

The authors are grateful to D. F. Agterberg, J. Akim- 
itsu, J. Flouquet, P.A. Frigeri, S. Fujimoto, J. Goryo, 
N. Hayashi, K. Izawa, N. Kimura, Y. Kitaoka, Y. Mat- 
suda, V. P. Minccv, H. Mukuda, E. Ohmichi, Y. Onuki, 
T. Shibauchi, R. Settai, T. Takeuchi, H. Tanaka, T. 
Tateiwa, H. Tou and M. Yogi for fruitful discussions. 
This study has been financially supported by the Nishina 
Memorial Foundation, Grants-in-Aid for Young Scien- 
tists (B) from the Ministry of Education, Culture, 
Sports, Science and Technology (MEXT), the Swiss Na- 
tionalfonds and the NCCR MaNEP. Numerical compu- 
tation was carried out at the Yukawa Institute Computer 
Facility. 

Appendix A: Linear Response Theory 

The dynamical spin susceptibility in the linear re- 
sponse regime is obtained by the Kubo formula as, 

-EE '[^*7(^ + 5"*' '')Sis{k, k + q)Gs{k + q)G^{k) 

7,<5 k 

-S>;,{k+ + q, k+)S';,{-k. - q, -k.)Fs{k + q)F^ik)] 

(A-l) 

where q = (q, iri„), k = [k, ia)„) and q is the momentum 
along the a6-plane. The spin operator S!^g{k + q, k) in the 
band basis has been given in eq. (17). 

Taking the limit f7„ and (f ^ 0, we obtain the 
uniform spin susceptibility Xnv = limg-_^olinif2„ 
which can be decomposed into a Van-Vleck and Pauli 
part as, 

X^fi. = limQ„^olim,-_^oX/:.i^(9), (A-2) 



We obtain the following expressions, 

7 k 

[S!^^{k, k)S;^{k, k)G^{k + q)G^{k) 
-Si;^{k+, k+)S;^{-k-,-k-)F^{k + 9)J;t(fc)], (A.4) 
Xl. = -Y.T.Wk,k)S^^{k,k)Gs{k)G,{k) 

-St;,{k+, k+)S;s{-k-,-k.)Fs{k)F^{k)]. (A.5) 

Assuming |A.y(fc)|, WF|<fH|/2 <C \a\ where vp is the 
Fermi velocity, the Van-Vleck part of spin susceptibility 
eq. (A.5) is obtained as in eq. (16). 

When we restrict to the AFM moment along the prin- 
cipal axis, namely Jiq \\ x,y or z, the relation U{—k) = 
e^^IU {k) holds with 6 an arbitrary phase factor. Here we 
denote 



/ = 



h 

±72 

1 
-1 



(A-6) 
(A-7) 



The sign of 72 in 7 is -|- for Kq \\ x, y and 
— for hq II z. According to eqs. (17) and (18), 
we obtain S!^^{—k,—k) ~ —St^^{k,k) for ^ = 

x,y and S^^{-k,-k) = S^^{k,k). If i;F|gH|/2 < 
jaj, the coefficient in eq. (A.4) is approximated as 
St;.^{k+, k+)S!;^{-k^,-k^) ^ St;^{k, k)S!^^{~k, -k) and 
the Pauli part of spin susceptibility is obtained as 
eq. (15). 

Appendix B: magnetic Properties in the d-wave 
State 

For the discussion for the intra-plane d-wave state we 
adopt the model eq. (29) but assume the tight binding 
parameters in eq. (27) as. 



{h,U,n) = (1,0.2,0.8), 



(B-1) 



P _ _ V 

^fj,!/ X/.it' 'X.fJ.l/' 



(A-3) 



with all other parameters zero. This parameter set leads 
to the nearly half-filled band with quasi-two dimensional 
Fermi surface and leads to the dominantly dx2_y2-wave 
SC state for the parameter set {A) U > Q,V = -Q.W. 
The order parameters arc described as $(fc) = 5 + 
T] cos fcx — cos ky with 6 = and rj = 1 at hq = h = 0. Our 
analysis confirms |^|, |1 — jyj <C 1. In general, the d^2_y2- 
wavc state is admixed with the fx{K^^y^)- and fy(^^2_y2y 
■ wave order parameters owing to the ASOC. However, the 
/-wave component does not appear in the mean field so- 
lution of the effective model eq. (29) because interactions 
beyond the nearest neighbor sites are neglected. 

We calculate the critical magnetic field /ic2 by solv- 
ing the linearized mean field equation eqs. (40) and (41) 
and show the result in Fig. 4. The spin susceptibility is 
calculated on the basis of eq. (43) by solving the mean 
field equation eqs. (35-38). In Fig. B.l we show the spin 
susceptibility below Tc for various magnetic fields. These 
results should be contrasted to those for the p-wave state 
(Figs. 4, 6, 10 and 11). 
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Fig. B-1. (Color online) The spin susceptibility along the [100]- 
and [010]-directions in the d-wave state without AFM order. We 
assume U > 0, V = -0.8U, J = 0, a = 0.3 and /iq = 0. The 
magnetic field is chosen as h = O.lTc, h = Tc, h = l.STc and 
h = 2Tc from the bottom to the top. 
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